DRAFT  SF  298 


1.  Report  Date  (dd-mm-yy) 

3.  Dates  covered  (from...  to  ) 

4.  Title  &  subtitle 

Control  of  Slowly  Varying  LPV  Systems:  An  Application  to 
Flight  Control 

5a.  Contract  or  Grant  # 

5b.  Program  Element  # 

6.  Author(s) 

5c.  Project  # 

5d.  Task  # 

5e.  Work  Unit# 

7.  Performing  Organization  Name  &  Address 

8.  Performing  Organization  Report  # 

9.  Sponsoring/Monitoring  Agency  Name  &  Address 

10.  Monitor  Acronym 

11.  Monitor  Report  # 

12.  Distribution/Availability  Statement 

Distribution  approved  for  Public  Release,  Distribution  Unlimited 

13.  Supplementary  Notes  TR  downloaded  from  a  WWW  unrestricted  URL. 

14.  Abstract 

19970121  161 

15.  Subject  Terms 

Security  Classification  of 

19.  Limitation 
of  Abstract 

Unlimited 

20.  #  of 
Pages 

21.  Responsible  Person 
(Name  and  Telephone  #) 

16.  Report 
Unclass 

17.  Abstract 
Unclass 

18.  This  Page 
Unclass 

DTIG  QUAItlTF  XlfSP: 


^OTED  \ 


CONTROL  OF  SLOWLY  VARYING  LPV  SYSTEMS: 
AN  APPLICATION  TO  FLIGHT  CONTROL  * 

Lawton  H.  Lee^  Mark  Spillman^ 

Dept,  of  Mechanical  Engineering  WL/FIGC-3,  Wright  Laboratory 

Univ.  of  California,  Berkeley  CA  94720  Wright-Patterson  AFB,  OH  45433-7531 


Abstract 

Recent  results  in  parameter- dependent  control  of  lin¬ 
ear  parameter- varying  (LPV)  systems  are  applied  to 
the  problem  of  designing  gain-scheduled  pitch  rate 
controllers  for  the  F-16  VISTA  (Variable- St  ability 
In-Flight  Simulator  Test  Aircraft).  These  meth¬ 
ods,  based  on  parameter-dependent  quadratic  Lya¬ 
punov  functions,  take  advantage  of  known  a  priori 
bounds  on  the  parameters’  rates  of  variation  (the 
bounds  may  themselves  be  parameter- varying) .  The 
controller  achieves  an  induced-£2-norm  performance 
objective;  Level  1  flying  qualities  are  predicted.  Sub- 
optimal  solutions  are  obtained  by  solving  a  con¬ 
vex  optimization  problem  described  by  linear  ma¬ 
trix  inequalities  (LMIs).  Incorporation  of  D-K  iter¬ 
ation  with  “constant  T)-scales”  provides  robustness 
to  time- varying  uncertainty.  Parameter- varying  per¬ 
formance  weights  are  used  to  shape  the  desired  per¬ 
formance  at  different  points  in  the  design  envelope. 


1  Introduction 

The  area  of  analysis  and  control  of  linear  parameter- 
varying  (LPV)  systems  has  received  much  recent 
attention,  primarily  in  order  to  develop  systematic 
techniques  for  gain-scheduling.  These  systems  re¬ 
semble  linear  systems  that  depend  on  one  or  more 
time- varying  parameters;  nonlinear  systems  are  of¬ 
ten  modelled  in  this  form  via  a  parameterized  family 
of  linearizations.  The  analysis  of  LPV  systems  dif¬ 
fer  from  that  of  linear  time-varying  (LTV)  systems  in 
that  it  considers  whole  families  of  parameter  trajec¬ 
tories;  moreover,  the  parameter  values  are  available 
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only  in  real  time,  not  in  advance. 

The  classical  approach  to  gain-scheduled  Tloo  con¬ 
trol  involves  designing  an  (LTI)  controller  for 
each  of  a  parameterized  family  of  linearizations  and 
then  interpolating  controller  gains  by  operating  con¬ 
dition.  This  heuristic  approach  yields  satisfactory 
results  if  the  parameters  are  sufficiently  “slowly- 
varying.”^^  Early  results  in  so-called  “LPV  synthe¬ 
sis”  explicitly  account  for  this  time- variation  using 
scaled  small-gain  arguments^’ or  single  quadratic 
Lyapunov  functions  (SQLF);^'®'®  those  designs  tend 
to  be  conservative,  though,  partly  because  they  al¬ 
low  the  parameters  to  vary  arbitrarily  quickly.  More 
recent  results^’ parameter-dependent  Lya¬ 
punov  functions  (PDLF)  to  factor  in  a  pnon  bounds 
on  the  parameters’  rates  of  variation,  reducing  this 
conservatism. 

In  this  paper  the  PDLF  technique^’ is  incor¬ 
porated  into  the  D-K  iteration  framework  in  order 
to  design  robust,  parameter- varying  pitch  rate  con¬ 
trollers  for  the  F-16  VISTA  (Variable  Stability  In- 
Flight  Simulator  Test  Aircraft) .  Sub-optimal  solu¬ 
tions  are  obtained  by  solving  a  convex  optimization 
problem  described  by  a  system  of  linear  matrix  in¬ 
equalities  (LMIs),  for  which  efficient  algorithms  are 
available.  Parameter- varying  performance  weights 
are  used  to  smoothly  vary  the  desired  performance 
within  the  design  envelope. 

The  remainder  of  this  paper  is  organized  as  follows. 
Section  2  reviews  some  results  on  the  control  of  LPV 
systems.  Section  3  provides  a  robustness  test  for 
LPV  systems.  Section  4  applies  these  techniques  to 
the  F-16  VISTA.  Section  5  summarizes  the  paper 
and  discusses  future  prospects. 

This  paper  uses  standard  notation.  In  addition, 
gnxn  (^gj3^otes  the  set  of  real,  symmetric,  n  x  n  ma¬ 
trices.  For  any  matrix  X  £  ^  X  >0  and  X  <  0 

respectively  denote  positive-definiteness  (all  of  its 
eigenvalues  are  positive)  and  negative-definiteness 
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(all  of  its  eigenvalues  are  negative).  Continuously 
differentiable  functions  are  called  and  C^[VyW) 
denotes  the  set  of  functions  from  V  to  W. 


2  Control  of  LPV  Systems 


This  section  reviews  some  results  on  the  control 
of  continuous-time,  linear  parameter- varying  sys¬ 
tems.  The  reader  may  consult  the  applicable  ref- 
erences*^’^^’^*^  (and  the  papers  cited  therein)  for  de¬ 
tails. 


An  LPV  system  G(s,  p)  is  a  finite-dimensional  linear 
system 


C{p{t))  D{p(t)) 


x{t) 

dit)  J 


(2.1) 


whose  state-space  data  are  known  continuous  func¬ 
tions  of  time- varying  parameters  denoted  by  />  E 
The  s  parameter  values  are  not  known  in  advance; 
rather,  they  are  measured  in  real-time.  Assume  that 
the  bounded  vector- valued  parameter  signal  p  E  Coo 
is  a  piecewise-C^  function  of  time  and  that  there  ex¬ 
ists  a  compact  set  'P  C  R^  for  which  p[i)  E  P  for  all 
t  >  0. 


Also  assume  that  the  rates  of  variation  of  the  first 
s  parameters  , . . . ,  ps  are  each  bounded  in  magni¬ 
tude  by  known  positive  scalars  . . . ,  i.e., 

\pi  {t)  I  <  z/,  for  alH  >  0  (z  =  1 , . . . ,  s) . 

Denote  these  s  rate-limited  parameters  and  rate 
bounds  by  p  :=  (pi,...,pj)  E  R^  and  v 
(z^i , , . . ,  z/j)  E  R^ ,  respectively.  Now  v  may  itself  be 
a  continuous  function  of  the  parameters  p,  so  that 
parameter  trajectories  must  obey  the  differential  in¬ 
clusion 


|Pi(<)|  <  n{p{i))  for  all  <>  0  (i  =  1, . . . ,  s) 

Parameter  trajectories  satisfying  the  above  condi¬ 
tions  for  given  P  and  v  will  be  called  allowable. 
Note  that,  at  the  cost  of  added  notation,  one  can 
expand  the  results  in  this  section  to  separate  upper 
and  lower  bounds  of  the  rates  of  variation. 


Lemma  2.2  Given  the  LPV  system  in  (2.1)  and  a 
performance  level  7  >  0,  suppose  there  exists  a  ma¬ 
trix  function  W  E  C^(R^,5”^”)  such  that  W[p)  >  0 
and  (omitting  dependence  on  p  and  p) 


A'^W  +  WA  +  Y:Li  w:  ^ 

B^W  '  -'fl 

C  D  -7/ 


<0 


(2.3) 

for  all  G  {—Vi{p),i'i{p)]  at  each  parameter  value 
p  E  P.  Then  for  any  allowable  parameter  trajec¬ 
tory  the  LTV  system  governed  by  (2.1)  is  exponen¬ 
tially  stable.  Furthermore,  there  exists  71  E  [0,7) 
for  which  ||e||2  <  7i||d||2  for  all  d  e  C2  (if  a)(0)  0). 


The  constraints  PV(p)  >  0  and  (2.3)  represent  con¬ 
vex  linear  matrix  inequality  (LMI)  constraints  on  the 
variables  W  and  7.  Although  these  LMIs  are  clearly 
infinite  dimensional,  one  can  compute  solutions  us¬ 
ing  the  approximate  method  described  in  the  sequel. 


2.2  Output-feedback  synthesis 

Consider  an  LPV  plant  in  the  standard  form 


X 

■  A{p)  B,{p)  B2{p)  ■ 

X 

e 

= 

Clip)  Diiip)  DM 

d 

-  y . 

C2{p)  D2i{p)  D22{p) 

u 

(2.4) 

where  x  E  R”,  e  E  R^®,  d  E  R”'^,  u  E  R”“,  and 
y  E  R”^,  and  other  quantities  are  dimensioned  ap¬ 
propriately.  By  assuming  regularity  (i.e.,  D12  and 
D21  are  full  rank)  and  Dn  =  0,  (2.4)  can  be  trans¬ 
formed  into^ 


~  X  ~ 

A(/?)  Pii(p)  Bi2{p)  B2{p) 

~  X  ’ 

ei 

Cn{p)  0  0  0 

di 

62 

Gi2(p)  0  0/ 

d2 

.  y  . 

C'2(p)  0  I  0  _ 

_  u  ^ 

(2.5) 

The  LPV  'y-Performance/u- Variation  problem  con¬ 
sists  of  finding  a  parameter- varying  controller 


Xc 

Ak{p,  p) 

Bk{p,  'p)  ' 

Xc 

u 

.  Ck(p,p) 

Bk{p,^)  . 

y 

(possibly  dependent  on  p)  for  which  the  closed-loop 
system  (omitting  dependence  on  p  and  p) 


2.1  Induced  £2“norm  Analysis 

Given  a  family  of  allowable  parameter  trajectories 
defined  by  a  parameter  set  P  and  a  rate-of- variation 
bound  z/,  one  can  bound  the  induced  C2  norm  of  an 
LPV  system  using  a  parameter-dependent  quadratic 
Lyapunov  function. 


X 

X 

Xc 

'o 

Xc 

ei 

G5 

di 

-  ^2  . 

_  ^2  _ 

satisfies  the  conditions  of  Lemma  2.2  for  a  desired 
closed-loop  norm  7  >  0.  The  following  theorem  gives 
necessary  and  sufficient  conditions  for  solvability. 
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Theorem  2.6  Given  the  compact  set  'P,  the  vector¬ 
valued  function  v{p),  the  scalar  7  >  0,  and  the 
open-loop  system  (2.5),  the  LPV  7- Perform ance/z/- 
Variation  problem  is  solvable  if  and  only  if  there  ex¬ 
ist  matrix  functions  X  G  and  Y  E 

such  that  (omitting  dependence  on  p 

and  p) 


B^iX 

L  Cl 

AY  YYA^ 

'  Z^i=i  dpi 
CiiY 

L  BJ 


X  I 
I  Y 

for  all  p  E  P,  where 

A  =  A  —  R12C2, 

A  —  A  —  R2C12, 


XBn 

Cf 

<  0 

-7/ 

0 

0 

(2.7) 

YCfi 

Bi 

<  0 

-yl 

0 

0 

-7^  . 

(2.8) 

>  0 

(2.9) 

Cf  —  [Cfi  C12] 

Bi  =  [Bn 

^12] 

The  notation  ili(-)i  indicates  that  the  inequalities 
must  hold  for  every  combination  of  +{')i  and  —(*)«. 
Therefore,  (2.7)  and  (2.8)  each  represent  2^  LMPs. 


Controllers  specified  by  arbitrary  parameter- varying 
solutions  X{p)  and  Y{p)  typically  depend  explicitly 
on  the  parameter  derivative  p.  On  the  other  hand, 
Becker^  shows  that  in  the  case  where  the  conditions 
of  Theorem  2.6  are  satisfied  when  X  or  Y  is  constant 
with  respect  to  p,  one  can  derive  formulas  for  strictly 
proper  {Dk  =  0)  controllers  that  are  independent  of 
p.  Specifically, 

Ak  =■  N  ^  -|-  X(^A  -|-  B2F  -j-  T(72)Y  T  Z'\h/[  ^ 

Bk  =  N-^XL 

Ck  =  FYM-'^ 

where 

F  =  -[jB^Y-^  +  Cu] 

L  =  -[jX-^C^  +  Bn] 

Z  =  [XiBiB'l +LBl^)  +  {ClC,  +  Cf^F)Y]/j 


and  M  &  AT  are  chosen  as  follows.  A  constant  X  and 
parameter- varying  Y  (p)  admit  the  controller  defined 
by  choosing 


M  =  I-  YX,  N  =  I 


Similarly,  a  parameter- varying  X(p)  and  constant  Y 
admit  the  controller  defined  by  choosing 

M  =  Y-^~X,  N  =  Y 

Holding  both  X  and  Y  constant  recovers  the  conser¬ 
vative  SQLF  controller,®’®  which  allows  for  arbitrar¬ 
ily  fast  parameter  variation. 


2.3  Computing  Solutions 


The  infinite-dimensionality  of  the  inequalities  in 
Theorem  (2.6)  demands  approximate  methods  of  so¬ 
lution  for  the  sake  of  practical  computation.  One 
such  method  follows:^®’ 

Pick  scalar  basis  functions  {/*  E  C^(R*,  R)}^^  and 
{dj  ^  and  search  over  those  X(p)’s 

and  Y(p)’s  that  are  linear  combinations 

x{p)  =  Er=i  fi  m, ,  Yip)  =  Efji  gj  im 


using  constant  matrices  {Xj  E  and  {Yj  E 

^nxnj.^^  Then  (2.7)-(2.9)  can  be  rewritten  as 
(omitting  dependence  in  p  and  p) 


ES  +  XiA)B  -  7C|'C2 

+  EU±-'=ESIb^. 

Cl 

E,7i  +  YA'r)gj  -  7S2BJ 


EjliSjCuYj 

BT 

'ES/.^.' 


E,tv,r,  J 


(*)  (^) 


-7/ 

0 


-7/ 

0 


-7I 


M  (^) 


-7J 


>  0 


<  0 

(2.10) 

<  0 

(2.11) 

(2.12) 


where  the  shorthand  (•^)  denotes  the  implied  trans¬ 
poses.  The  problem  thus  consists  of  finding  real, 
symmetric  matrices  {XJ^^  and  {Yj}^^  that  sat¬ 
isfy  the  above  inequalities  for  all  p  E  P,  still  an 
infinite-dimensional  problem. 


Now  approximate  the  parameter  set  'P  by  a  grid  of 
L  points  {pk  E  defining  {pk  E  R^}^=i  ac¬ 

cordingly,  and  solve  the  inequalities  (2.10)-(2.12)  at 
these  grid  points.  The  conditions  represent  up  to 
1/(2®+^  +  1)  LMIs  in  the  Nx  +  Ny  matrix  variables 
{{Xf}^^,  Note  that  minimizing  7  is  a 

convex  optimization  problem,  since  the  inequalities 
are  affine  in  7  as  well. 
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3  Robust  LPV  Systems 


4 


A  Design  Example 


This  section  derives  a  simple  robustness  test  and 
proposes  an  ad  hoc  method  for  making  the  preced¬ 
ing  control  synthesis  robust.  Consider  the  analysis 
of  an  LPV  system 


put  in  feedback  with  a  block- diagonal,  memoryless 
linear  time-varying  (uncertainty)  operator 

w{t)  =  A(^)2:(t) 


This  section  presents  an  application  of  LPV  syn¬ 
thesis  to  the  pitch  rate  control  of  the  F-16  VISTA 
over  a  specified  flight  envelope;  a  similar  problem 
has  been  addressed^^  using  a  small-gain  method. 
The  physical  plant  and  certain  design  weights  are 
parameter- varying,  and  a  D-K  iteration-like  process 
is  employed  to  enhance  robust  performance. 

4.1  Plant  Modelling 

This  design  uses  the  standard  short-period  equations 
of  motion^ 


where  w,z  G  d  E  and  e  E  R^®.  Assume 
that  the  block  structure  of  A  is  consistent  with  a  set 
A  C  Define  a  corresponding  set  of  scalings 


Sa  :=  {S  E  :AS  =  SA  for  all  A  E  A} 


An  elementary  small-gain  argument  establishes  the 
following  sufficient  conditions  for  robustness  to  (ar¬ 
bitrarily  quickly)  time- varying  uncertainty. 

Proposition  3.1  Suppose  there  exists  a  continuous 
matrix  function  5  :  R^  Sa  such  that  for  any  al¬ 
lowable  parameter  trajectory  p  :  R  -4  R®  the  result¬ 
ing  LTV  system 


Gs{s,p) 


S{p)  0 
0  In. 


G{s,p) 


s-Hp)  0 

0  ^rid 


is  exponentially  stable  and  ||G5(s5  p)||j2  <  7  (given 
zero  initial  conditions) .  Then  for  any  allowable  pa¬ 
rameter  trajectory  p  and  any  A  satisfying 


a(A(f))  <  1/7 


the  LTV  system  Fu{G{s^  p).  A)  also  is  exponentially 
stable  and  ||Fti(G(s,  p),  A)||j2  <  7  for  zero  initial 
conditions. 


This  suggests  that  one  can  design  robust  LPV  con¬ 
trollers  by  alternating  (in  obvious  fashion)  between 
controller  synthesis  and  computation  of  parameter- 
varying  scaling  matrices  5(/>),  as  in  D-K  iter¬ 
ation.  The  design  example  in  this  paper  uses 
an  ad  hoc  choice  of  S[p):  “zeroth-order”  fits  to 
frequency-dependent  D-scales  that  are  obtained  at 
each  “frozen”  grid  point  by  applying  the  appropriate 
p-Tooh  commands  for  closed-loop  analysis.  There  is 
also  a  more  rigorous  iterative  approach^  that  will  not 
be  addressed  in  this  paper. 


a 

Za{p)  1 

.  q  . 

_  M„(p)  M,{p) 

P  =  P={M,h) 


Zk{p) 

MsAp) 


Se 

(4.1) 


ignoring  the  aerodynamic  effects  of  the  trailing  edge 
flaps.  Only  the  longitudinal  dynamics  of  the  aircraft 
are  considered;  the  roll,  yaw,  and  sideslip  angles  are 
assumed  to  be  zero.  In  (4.1)  the  states  (a,g),  input 
Se ,  and  parameters  (M,  h)  respectively  denote  angle- 
of- attack  &  pitch  rate,  elevator  deflection,  and  Mach 
number  &  altitude. 


The  Flight  Dynamics  Directorate  of  the  Wright  Lab¬ 
oratory  uses  a  high-fidelity,  six  degree-of-freedom, 
nonlinear  model  to  simulate  the  F-16  VISTA. ^  This 
simulation  model  includes  accurate  descriptions  of 
the  propulsion  system,  actuators,  sensors,  distur¬ 
bances,  payload,  atmosphere,  rigid-body  equations 
of  motion,  and  aerodynamics  for  a  wide  range  of 
Mach  numbers,  altitudes,  and  angles  of  attack. 

The  LPV  short-period  model’s  state  space  data  (the 
dimensional  coefiScients  Za^  Mq,  Zs^,  and  Ms^) 
were  obtained  by  trimming  and  linearizing  the  non¬ 
linear  model  at  level  flight  for  the  flight  conditions 
marked  in  Table  1.  The  design  region  V  was  chosen 
accordingly,  and  these  21  grid  points  were  used  for 
the  controller  synthesis. 


h\M 

0.35 

0.45 

0.55 

0.65 

0.75 

0.85 

25000  ft 

X 

X 

X 

X 

15000  ft 

X 

X 

X 

X 

X 

5000  ft 

X 

X 

X 

X 

X 

X 

1000  ft 

X 

X 

X 

X 

X 

X 

Table  1:  Grid  points  used  for  modeling  &  synthesis 


Data  on  the  excess  thrust  and  rate-of-climb  of  the 
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F-16  VISTA  suggest  the  bounds 


Mp)  =  2(i%l^  +  M2|«|)  >  \M{t)\ 
V2{p)  =  a{h)M  >  |A(<)| 

on  the  parameters’  rates  of  variation,  where  a{h) 
denotes  the  speed  of  sound  as  a  function  of  altitude 
and  Vt  =  a[h)M  denotes  the  aircraft’s  true  velocity. 
Note  that  these  bounds  are  conservative;  achieving 
them  would  require  vertical  flight,  for  example. 


4,2  Problem  Setup 


The  objective  here  is  to  design  for  the  F-16  VISTA  a 
pitch-rate  controller  that  provides  robust  command 
tracking  with  predicted  Level  1  handling  qualities. 
Time-domain  specifications  for  pitch-rate  response 
are  illustrated  in  Fig.  1  and  listed  in  Table  2.  Note 
that  the  ^^rise-time”  parameter  At  varies  with  the 
true  velocity  Vt  (in  ft/sec). 


meets  Level  1  flying  qualities  over  the  entire  design 
envelope.  The  first-order,  parameter- varying  com¬ 
mand  and  performance  weights 


Wr{s,p) 

W,{s,p) 


s  +  100  10 

100  5+10^”'''^^^^ 

^  +  80  4  1 

80  s  + 4  0.05gmaa7(p) 


reflect  a  uniform  10  rad/s  command  bandwidth, 
a  maximum  pitch-rate  command  qmax  that  varies 
across  the  design  envelope  as  shown  in  Table  3,  and 
steady-state  tracking  within  5%. 


h\M 

0.35 

0.45 

0.55 

0.65 

0.75 

0.85 

25000  ft 

13 

17 

20 

22 

15000  ft 

12 

15 

18 

21 

23 

5000  ft 

10 

13 

16 

19 

22 

24 

1000  ft 

11 

14 

17 

20 

23 

25 

Table  3:  Max.  pitch-rate  command  Qmax  (deg/s) 


The  design  plant  also  includes  actuator  dynamics, 
additive  sensor  noise,  multiplicative  input  uncer¬ 
tainty,  and  penalties  on  the  elevator  deflection  angle 
and  rate.  The  first-order  actuator  model 


Ga{s) 


20.2 
s  20.2 


reflects  a  bandwidth  of  20.2  rad/s. 
weight 


21  0 
0  70 


The  control 


reflects  elevator  deflection  angle  and  rate  limits  of 
21  deg  and  70  deg/s,  respectively.  The  parameter- 
dependent  noise  weight 


Figure  1:  Pitch  rate  handling  qualities  specifications 


Parameter 

Level  1 

Level  2 

max  it  1  (sec) 
max  Aq2/Aqi 
max  At  (sec) 
min  (sec) 

0.12 

0.30 

m/Vx 

9.0/Vt 

0.17 

0.60 

1600/7t 

3.2/I^t 

Table  2:  Pitch  rate  handling  qualities  specifications 


Wn(p)  = 


0.5  0 

0  O-O^qmaxip) 


anticipates  0.5  deg  of  measurement  error  for  a  and 
0.03  grnarc  deg/s  (about  3%)  for  g.  The  uncertainty 
weight 

Wu  =  0.1 


represents  10%  parametric  and/or  dynamic  model¬ 
ing  error. 


This  design  uses  the  model-matching  control  struc¬ 
ture  shown  in  Fig.  2.  The  second-order  reference 
model 


Gref  (5) 


u  =  4  rad/s,  C  =  0.6,  1/T  =  10  rad/s 


4.3  Design  Results 


5 


Several  D-K  iterations  (using  LPV  synthesis)  were 
performed  in  MATLAB  on  the  7th-order  design 
plant.  LMILab^o  was  used  to  solve  the  controller 


Figure  2:  The  design  plant 


synthesis  LMIs,  and  /i-Tools^  was  used  for  closed- 
loop  robustness  analysis.  The  elementary  basis  func¬ 
tions 

[flip)  f2ip)  fsip)]  =  [giip)  92{p)  gziP)]  =  [lMh] 

selected  according  to  previous  experience,^’ 
were  used  to  vary  the  matrices  X{p)  and  Y{p)  with 
various  degrees  of  complexity:  constant  X  h  Y y 
varying  X[p)  &  constant  Y,  constant  X  &  varying 
Y{p)y  and  varying  X{p)  kY{p). 


D-K 
Iter.  # 

7 

X,Y 

7oo 

7 

X,Yip) 

Too 

1 

2.60 

1.90 

1.79 

1.49 

1.14 

1.13 

2 

2.43 

1.75 

1.73 

1.22 

1.05 

1.04 

3 

2.42 

1.74 

1.72 

1.19 

1.04 

1.03 

4 

2.42 

1.74 

1.72 

1.18 

1.04 

1.03 

Table  4:  Closed-loop  performance  levels 


The  closed-loop  time-varying  performance  level  7 
achieved  using  Theorem  2.6  is  shown  in  Table  4;  con¬ 
straining  the  rate  of  variation  of  M  an  clearly  of¬ 
fers  a  significant  improvement  in  performance.  Also 
included  are  the  maximum  “frozen- point”  (i.e.  gain- 
scheduled)  norms 

7oo  =  m^||Gc/p(s,/?)||oo 

p€T 

and  structured  singular  values  (cf.  Figure  3) 
p  =  m^  sup  p[Gcip{3U)y  p)] 

obtained  via  pointwise  Woo  and  //-analysis  of  the 
closed-loop  systems.  These  indicate  the  controller’s 
robustness  at  constant  flight  conditions.  The  two 
designs  with  varying  X{p)  offer  no  appreciable  im¬ 
provement  over  the  corresponding  constant-X  de¬ 
signs,  so  they  have  been  omitted. 


Figure  3:  Closed-loop  p{ju})  after  4th  iteration 

4.4  Nonlinear  Simulation 

Figure  4  shows  the  results  of  a  high-fidelity, 
parameter-varying,  nonlinear  simulation  of  the 
closed-loop  step  response,  which  demonstrates  pre¬ 
dicted  Level  1  flying  qualities.  The  aircraft  is  ini¬ 
tially  perturbed  from  trimmed,  level  flight  at  M  = 
0.75  and  h  =  5000  ft.  The  response  of  the  constant- 
Xy  Y  controller  is  also  shown,  for  comparison. 


5  Conclusions 


Recent  results  in  parameter-dependent  control  of  lin¬ 
ear  parameter-varying  systems  are  applied  to  the 
problem  of  robust,  gain-scheduled  pitch  rate  con¬ 
trol  for  the  F-16  VISTA.  Using  parameter-dependent 
Lyapunov  functions,  a  priori  bounds  on  the  pa¬ 
rameters’  rates  of  variation,  LMI-based  convex  op¬ 
timization,  and  parameter- varying  design  weights, 
this  method  achieves  an  induced-£2“norm  perfor¬ 
mance  objective  while  predicting  Level  1  flying  qual- 
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'igure  4:  Parameter- varying  nonlinear  simulation 


ities  throughout  the  design  envelope.  Straightfor¬ 
ward  D-K  iteration  with  “constant  £)-scales”  pro¬ 
vides  robustness  to  time-varying  uncertainty.  On¬ 
going  research  includes  expanding  the  parameter  set 
(to  encompass  the  full  flight  envelope)  and  including 
a  parameter-varying  reference  model  (while  main¬ 
taining  Level  1  flying  qualities),  and  incorporating 
the  robustness  scales  into  the  synthesis  LMIs. 
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